Model dependence of single-energy fits to pion photoproduction data 
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Model dependence of multipole analysis has been explored through energy-dependent and single- 
energy fits to pion photoproduction data. The MAID energy-dependent solution has been used as 
input for an event generator producing realistic pseudo data. These were fitted using the SAID 
parametrization approach to determine single-energy and energy-dependent solutions over a range 
of lab photon energies from 200 to 1200 MeV. The resulting solutions were found to be consistent 
with the input amplitudes from MAID. Fits with a x-squared per datum of unity or less were 
generally achieved. We discuss energy regions where consistent results are expected, and explore 
the sensitivity of fits to the number of included single- and double-polarization observables. The 
influence of Watson's theorem is examined in detail. 
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I. INTRODUCTION 

The non-perturbative regime of quantum chromody- 
namics (QCD) is characterized by an array of hadronic 
resonances. The gross features of this spectrum have 
been understood theoretically within the context of con- 
stituent quark models. Experimentally, measurements of 
reaction observables from an array of collision processes, 
in particular pion photoproduction, have shown the limi- 
tations of the quark model description. Precision electro- 
magnetic facilities around the world have begun to give 
a more detailed picture of the hadronic resonances, par- 
ticularly for the resonances of the nucleon, and ushered 
in a renaissance in the field of hadronic reaction theory. 

Crucial to this program is the determination of reac- 
tion amplitudes from experimental observables. Model 
independent extraction of reaction amplitudes in both 
spin (eg. CGLN, helicity, transversity, ...) and partial- 
wave bases is a well-studied yet complex task. Typi- 
cally, coupled-channel models based on Lagrangians of 
hadronic effective field theory involve hundreds of bare 
parameters in order to obtain realistic descriptions of 
data covering the first, second, and third resonance re- 
gions. Complementary to this approach is the extraction 
of amplitudes with smaller numbers of parameters and 
minimal model dependence. 

Two types of analysis, typically performed in 
parametrizing the reaction amplitudes, lead to energy 
dependent (ED) or global solutions and energy indepen- 
dent or single-energy (SE) solutions. The SE amplitudes, 
generally extracted through the partial-wave analysis of 
scattering or reaction data, are often used as the start- 
ing point for more involved multi-channel analyses [T]. 
The term data is commonly applied to these sets of am- 
plitudes, implying a relatively model-independent link to 
the underlying experimental data. This approach has 



been successful for some hadronic and nuclear collision 
processes. For example, complete or nearly complete sets 
of data have been assembled over restricted kinematic re- 
gions for nucleon-nucleon [2] and pion-nucleon [3| elastic 
scattering. These have allowed a formal study of the 
ambiguities associated with the amplitude reconstruc- 
tion process and yielded a set of useful amplitudes [4] em- 
ployed, for example, in the determination of the isospin 
\ N and isospin- 1 A nucleon resonances [5]. 

In pion photoproduction, which is the focus of the 
present study, SE solutions are generally determined from 
data within a narrow range of energies, with less than 
the number of observables required for a direct, unique 
amplitude reconstruction. In order to obtain a stable so- 
lution, further constraints are required. Generally, the 
SE solutions or fits are constrained by the results of a 
global ED fit. Below, we will discuss how the ED and 
SE fits are related, and how this procedure differs from 
amplitude reconstruction. 

Multipole fits from the MAID and SAID groups have 
been updated to include improved measurements as they 
have become available. Both ED and SE results have 
been presented and some multipoles show significant dif- 
ferences. In order to avoid the influence of systematic un- 
certainties in the database, as a potential source of these 
discrepancies, we have generated pseudo data based on 
the MAID and SAID ED fits. Multipoles extracted from 
these idealized datasets show much less variation than is 
seen in fits to the existing experimental datasets. Below 
we compare multipoles analyzed using the two different 
techniques, ED and SE, and suggest energy ranges over 
which varying degrees of model independence can be ex- 
pected. 

In the next section, we provide information sufficient 
to fix the phase and sign conventions of amplitudes and 
observables. It is important to note that the phase and 
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FIG. 1: Frames for polarization vectors in the CM. 



sign conventions adopted by the MAID and SAID groups 
have been consistent for many years and, as such, have 
been adopted by a large number of experimental groups. 
Recent publications [HI [Zj, however, have departed from 
this convention and we offer the complete set of equations 
in Sec|n]in order to facilitate comparisons. 

Section |rPT| dcscribes the approaches taken by the SAID 
and MAID groups in the determination of SE solutions. 
In the following section, Sec. |IV[ we describe the event 
generator used to produce the pseudo data employed in 
the SE and ED fits. The results of our fits are given 
in a comparison of observables and multipole amplitudes 
in Sec(V] Our ED and SE solutions span a range of lab 
photon energies E 1 from 200 to 1200 MeV or total CM 
energies W from 1120 to 1770 MeV. In order to limit 
the length of this article, we have included figures for 
lab photon energies of 340 and 600 MeV corresponding 
to CM energies of 1230 and 1420 MeV. The lower en- 
ergy is in the region dominated by the elastic P33 partial 
wave. The higher energy is characterized generally by 
the onset of larger inelasticities, except in the P33 partial 



to the exclusive photoproduction cross section 



wave. Section VI gives a discussion of the results and our 
conclusions. 



q f M N 
k \4irW 



(2) 



where q{k) is the magnitude of the three-momentum of 
the pion (photon), is the mass of the nucleon, W 
is the CM scattering energy, and t^L is the transition 
matrix element for pion photoproduction; the spin and 
internal quantum numbers in this expression are implied. 
This matrix element is simply expressed in terms of the 
Pauli spinors of the initial and final state nucleons. 



A. CGLN and Helicity Amplitudes 



II. AMPLITUDES AND OBSERVABLES 



Expressed in terms of two-dimensional Pauli spinors, 
the matrix element of the electromagnetic current takes 
the form 



In general, pion photoproduction from the nucleon, 
7A — > 7rA, is expressed by four invariant amplitudes, 
Ai, which are covariant functions of two kinematical vari- 
ables, such as the Mandelstam variables s,t, or more of- 
ten E 7 , 8, the photon laboratory energy and the pion CM 
angle. For given initial and final spin states i, f the tran- 
sition matrix element is given by 



k=l 



(1) 



fi 4irW t 



(3) 



where Xf(i) i s a Pauli spinor for the nucleon in final 
(initial) state. Following Ref. [8 , the operator T is de- 
composed into four component amplitudes, Fi, called the 
"CGLN" amplitudes, as 



where Ufu-\[p) is the Dirac spinor of the final (initial) 
nucleon, normalized as u(p)u(p) — 2Mn, e M is the polar- 



ization vector of the photon, and Mjr are a set of Dirac 
matrices defined to be consistent with Ref. 8 . This leads 



-f 7 M 



= i(a-e)F 1 + (a-q)(axk)-eF 2 
+ i (e ■ q) (a ■ k)F 3 + i(e ■ q)[ff ■ q) F 4 



(4) 
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FIG. 2: Pseudo data (black circles) compared to published experimental data (open circles) and the MAID ED solution MD07 
(solid curve) at photon beam energies of 320 MeV for da/dD. (a), E (b) and T (c) and 340 MeV for the double spin observable 
G (d). 



where e M = (0, e) and e- k = for real photons. The mul- 
tipole series of the CGLN amplitudes takes the form[51[H]: 



F 1 =Y j [(lM l+ + E l+ )P l+x {x) 

1=0 

+{{l + l)M l _+E l _)P l _ 1 {x) 

OO 

F 9 =Y i [(! + l)M l+ +lM l _]Pl(x), 



(5) 
(6) 



1 = 1 
OO 



F 3 = j2m + -M l+ )p;' +1 ( X ) 
i=i 

+(Ei-+M l -)Pi- 1 (x)}, (7) 

OO 

F 4 = J2 Wi+ - E l+ - Ml -E^} Pi (x) , (8) 



1=2 



where x = cos 9 is the cosine of the scattering angle. 

The representation of the photoproduction amplitudes 
given by Eqs.([4])-(|8]) is useful in determining a consistent 
notation since the CGLN amplitudes take this form in 
the literature universally, to our knowledge. This pro- 
vides a context for the discussion of sign or phase con- 
ventions and various linear combinations of amplitudes 
to form different bases. In this way, one may readily dis- 
criminate various sign conventions used for observables. 



While more compact expressions for the observables of 
Eqs.(17|-(32), can be given in terms of the helicity or 



transversity amplitudes (see Table |TJ) , these amplitudes 
are not uniquely defined. Differing phases and varying 
linear combinations to define the helicity and transver- 
sity bases have been used in the literature. In this work, 
we use the conventions adopted by Walker [TO] for the 
helicity amplitudes. These conventions are also used in 
Refs.[TTJ[r2] of Barker, Donnachie, Storrow (BDS). The 
relations in Table [I] give the reaction observables explic- 
itly in terms of the helicity amplitudes using the conven- 
tions of Ref. [TU] . Comparisons are made there to works 
using differing sign conventions for convenience. 

The helicity amplitudes are related to the CGLN am- 
plitudes as 



Hi = — - f=. sin 9 cos - (F 3 + F 4 ) 
v 2 2 



(9) 



H 2 = V2 cos °- [(F 2 -F 1 )+ l C 2 ° Sd (F 3 - F 4 )] , (10) 



H 3 = sin 9 sin - {F 3 - F 4 ) , 



(11) 



H 4 = V2 sin °- \{F X + F a ) + ^|°^(F 3 + F 4 )] , (12) 

which are identical to Eqs.(24) of Ref.[in|. The outgoing 
pion direction is specified by the scattering angle 9 and 
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FIG. 3: The SE4 fit (solid) and ED4 fit (dashed) compared to pseudo data (open circles) generated from MD07 for the 7r°p 
channel for observables da/dD. (a), P (b), E (c), and T (d), at E 1 — 340 MeV. Each figure includes two, usually overlapping, 
curves. 



azimuthal angle </>. We have chosen 4> = by aligning the 
coordinate system with the reaction plane as in Fig.(fl]) 
(see the discussion immediately following Eq.( 16 ) below). 
The same helicity amplitudes are employed by BDS |lll 
I12j with, however, a distinct naming convention: 

S! = H lt N = H 2l 

D = H 3 , S 2 = (13) 

B. Coordinate Frames and Polarizations 

Experiments with three types of polarization may be 
performed in meson photoproduction: photon beam po- 
larization, polarization of the target nucleon, and polar- 
ization of the recoil nucleon. The target polarization is 
described in the frame {x : y,z} of Fig. ([!]), with the z- 

axis coincident with the photon momentum fc, the y-axis 
perpendicular to the reaction plane, y — k x qj sin0, and 
the x-axis given by x = y x z. For recoil polarization 
we use the frame {x' ,y' , z'}, with the z'-axis coincident 
with the momentum of the outgoing meson q, the y'-axis 
coincident with the y-axis, as for the target polarization, 
and the x'-axis given by x' = y x z* . 

The photon polarization is either linear or circular. 
A linearly polarized photon is defined by the degree of 



transverse polarization, Pt, and the angle ip of the po- 
larization plane relative to the reaction plane or, cquiv- 
alcntly, to the x direction. For example, a beam com- 
pletely polarized in the reaction plane has ip = (with 
Pt = 1); for a beam polarized perpendicular to the re- 
action plane, y, the polarization angle is ip = n/2 (with 
Pt = 1). For a photon of right- (left-) circular polarization 
P Q = ±1. 

The target nucleon polarization is specified by three 
polarization components (P x ,P y ,P z ) with respect to 
the (x,y,z) coordinate system, displayed in Fig. ([I]). 
The recoil nucleon polarization is specified similarly by 
(P x i , P y i , P z i) with respect to (x', y*, z") also shown in 
Fig.0. 



C. Polarization Observables 

We may classify the differential cross section for gen- 
eral polarization states of the beam, target, and recoil 
particles by three classes of double polarization exper- 
iments. Using the notation described in the previous 
subsection, measurements with polarized photons and a 
polarized target, the beam-target experiments, are given 
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FIG. 4: The SE4 fit (solid) and ED4 (dashed) compared to pseudo data (open circles) for n°p at 600 MeV, as in Fig.|3|. 



by the differential cross section 

1 - P T £cos2<p 



1 da 
(To dfl 



P x {-P T Hsm2p + P Q F) 
P y (T- P T P cos 2ip) 
P z (P T Gsm2^~P E); 



(14) 



for polarized photons and measured recoil polarization, 
beam-recoil experiments, we obtain 



1 da 
a dfl 



= 1 - P T T,cos2ip 

+ P x . {-P T O x . sin2<^ - P Q C X .) 

+ P yl (P — PtT cos 2ip) 

+ P z ,(-P T O z ,sm2i P -P Q C z ,)- 1 (15) 



for the polarized target and recoil polarization measure- 
ments, target-recoil experiments: 



1 da 



l + P y T + P y ,P + P x , {P X T X , - P Z L X ,) 
P y ,P y Z + P z ,(P x T z ,+P z L z ,). (16) 



In these equations, ag denotes the unpolarized differen- 
tial cross section. Here ip is the azimuthal angle of the 
photon polarization vector with respect to the reaction 
plane. Alternatively, one could fix the photon polariza- 
tion vector and observe the pion out of the polarization 
plane of the photon. In this case one would obtain a 
pion angle <fi = —ip. This leads to a minus sign in all 
terms proportional to sin2(p. Equations (14 1, (15), and 
^ are identical with Eqs.(2), (3), and (4) of Ref.Ql] 
upon conversion of their ai — > Pi> , Oi — > Oi> , Ci — > CV , 
Ti — > Tii, and Li — > Li*. Note, however, that our coordi- 
nate frames are identical to those of BDS. The triple po- 
larization cross section, which specifies the polarizations 
of the beam, target, and recoil nucleon simultaneously, 
provides no additional information in the pion photopro- 
duction reaction. 



The spin observables, expressed in terms of the CGLN 
amplitudes, within the sign convention of Barker, Don- 
nachie and Storrow, have been consistently employed in 
both the MAID and SAID parametrizations. They are 
given explicitly as: 
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FIG. 5: The SE4 fit (solid) and ED4 (dashed) compared to pseudo data (open circles) for 7r°p beam-target observables, E (a), 
F (b), G (c), and H (d) at 600 MeV. 



t 
f 
P 
E 
F 
G 
H 
G x i 

c z , 
6 X , 
d z , 

L x ' 
L x > 



Re {F^Fx + F*F 2 + sin 2 9 (F 3 *F 3 /2 + F 4 *F 4 /2 + F 2 *F 3 + F*F 4 + cos 9 F 3 *F 4 ) - 2 cos 9 F*F 2 } p (17) 

- sin 2 9 Re {{F*F 3 + F 4 *F 4 ) /2 + F 2 *F 3 + F*F 4 + cos F 3 F 4 } p (18) 
sin 6> Im {F*F 3 - F*F 4 + cos 9 (F*F 4 - F 2 *F 3 ) - sin 2 9 F 3 *F 4 } p (19) 

- sin 9 Im { 2F X * F 2 + F* F 3 - F 2 * F 4 - cos 9 (F 2 * F 3 - F* F 4 ) - sin 2 F 3 * F 4 } p (20) 
Re {F^Fi + F*F 2 - 2cos0F 1 *F 2 + sin 2 9 (F 2 *F 3 + F*F 4 )} p (21) 
sin 9 Re {F* F 3 - F 2 F 4 - cos (F 2 * F 3 - Fj" F 4 ) } p (22) 
sin 2 6» Im {F*F 3 + F^FJ p (23) 
sin 9 Im {2F!*F 2 + F X *F 3 - F 2 *F 4 + cos 9 (F?F 4 - F 2 *F 3 )} p (24) 
sin 6> Re {F* Fi - F* F 2 - F* F 3 + F* F 4 - cos (F 2 * F 4 - Fj" F 3 ) } p (25) 
Re { 2F* F 2 - cos (F* F 1 + F* F 2 ) + sin 2 9 (F* F 3 + F 2 * F 4 ) } p (26) 
sin Im { F 2 * F 3 - Fj* F 4 + cos 9 (F 2 * F 4 - FiF 3 )} p (27) 

- sin 2 9 Im { F* F 3 + F 2 * F 4 } p (28) 

- sin Rc {F*Fi - F*F 2 - F 2 *F 3 + F*F 4 + sin 2 (F 4 *F 4 - F 3 *F 3 )/2 + cos 9 (F*F 3 - F 2 *F 4 ) } p (29) 
Re {2F 1 *F 2 - cos 9 {F*F 1 + F*F 2 ) + sin 2 (F*F 3 + F 2 *F 4 + F*F 4 ) + cos sin 2 9 (F*F 3 + F 4 *F 4 )/2} p (30) 

- sin 2 Rc {F*F 3 + F*F 4 + F 3 *F 4 + cos (F 3 *F 3 + F 4 *F 4 )/2} p (31) 
sin 9 Re {F 1 *F 4 - F 2 *F 3 + cos 9 (F?F 3 - F 2 *F 4 ) + sin 2 9 (F 4 *F 4 - F 3 *F 3 )/2} p, (32) 

I 



where £ = Sctq, etc. and p = q/k. 



III. SINGLE-ENERGY SOLUTIONS 



The SE solutions afford the most faithful description 
of the data, in terms of the lowest possible x-sq uare d, 
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TABLE I: Spin observables expressed by helicity amplitudes in the notation of Walker 10 . The sign definition is taken from 
Barker, Donnachie and Storrow [12] by replacing JV — > H 2 , Si — > Hi, S2 — > H4, D — > H3 as given in the text. This sign 
definition is used by SAID and MAID. In column three, we compare the sign definitions of Fasano, Tabakin and Saghai [13] . 
which is also used recently by Refs. [7] [14]. Furthermore we also give the sign definition of Chiang and Tabakin |15) and the 
relations to the Q observables, that are the basis of the Fierz consistency relations given in Ref. |15j . 



consistent with constraints from the ED solutions. They 
have been used to probe the data for structure that may 
not have been properly encoded by the smooth ED solu- 
tion. Both MAID and SAID SE solutions are determined 
by constrained fits to the data in sufficiently narrow en- 
ergy bins. These constraints limit the variability in the 
SE results through the assumption that the global ED fit 
is close to the 'true' solution. This assumption appears 
to be validated by the fact that the %-squared of both 
the ED and SE solutions suggest a realistic description 
of the data. Since the ED results are used to constrain 
the SE solutions in different ways within the MAID and 
SAID parametrization approaches, we begin this section 
with a brief description of their respective methods. 



A. SAID SE Solutions 

The SAID SE solutions are obtained in bins of the scat- 
tering energy spanning a few MeV. Starting values for 
the SE solutions are given by the multipole moduli and 
phases determined by ED solutions. In order to account 
for the variation of the modulus and phase over the width 
of an energy bin, a linear approximation to the energy de- 
pendence for each quantity is taken from the ED fit. The 
bin width increases from 5 to 20 MeV as the energy in- 
creases into regions with sparse measurements. Only the 
central values of the modulus and phase are searched for 
each SE solution; the slopes are held fixed. For the pur- 
poses of the present study, however, the pseudo data are 
generated at single energies coinciding with the central 
energy value of the bin. This removes the need to esti- 



mate the energy variation based on an ED fit. In order 
to constrain the SE multipole solutions to be close to the 
ED values, the ED multipoles themselves are ascribed 
arbitrary errors and fitted along with experimentally de- 
termined reaction data of the bin under consideration. In 
general, these multipole 'pseudo data' contribute a neg- 
ligible amount to the resultant chi-squared. 

In all previously published sets of SE fits, the multipole 
phases were fixed to ED values when fitting to the data. 
As a result, multipole phases would change only through 
a modification of the ED fit. The number of searched 
multipoles has generally been subjective but may be in- 
creased until a % 2 /data near unity is approached. 

As previously mentioned, the SAID SE fits were gen- 
erated originally to search for structures not seen in the 
ED fit, and they have been used to re-initialize the ED 
fit, with the hope that this procedure would converge to 
the correct solution. As described above, they are not 
independent of the ED fit. It should be emphasized that 
the ED and SE solutions are fits to the experimentally 
observed reaction data; the ED solution is not a fit to 
the SE values. 



B. MAID SE solutions 

The MAID partial wave analysis follows a similar two- 
step approach [17]. In the first step, a global ED solu- 
tion is determined first by fitting all experimentally ob- 
served reaction data, in a similar fashion to the SAID 
ED approach, in the range 140 MeV< E 1 <1610 MeV. 
This allows the determination of the phases of each mul- 
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tipole, i.e., the ratio Imt% /Ret% , where a is the mul- 
tipole, above the two-pion threshold. At energies be- 
low the two-pion threshold, this phase is constrained, 
by Watson's theorem, to be equal to the irN scatter- 
ing phase. In the second step we perform local SE fits 
to the data in energy bins of 10 MeV (over the range 
140 MeV< E 1 <460 MeV) and 20 MeV (for the higher 
energies) by varying the absolute values of the multipoles 
but keeping the phase fixed. In order to damp strong lo- 
cal variations we introduce a penalty factor similar to a 
Bayesian approach and minimize the modified x 2 func- 
tion 



N da 



E 



6j-e 



2 n W Vx< _jNa 



3 



A 



(33) 



The first term on the right-hand side of this equation is 
the standard x 2 function with 0^ the calculated and Q^ xp 
the measured observables, SQi the statistical errors, and 
Ndata the number of data points. In the second term, 
Nmuit is the number of the varied multipoles and Xj is 
the fitting parameter describing the deviation from the 
global fit, that is, the fitting procedure starts with the ini- 
tial value Xj = 1 corresponding to the global solution, for 
each multipole indexed by j. The quantity A enforces a 
smooth energy dependence of the single energy solution. 
In the limit of A — > oo we obtain the standard x-squared, 
x1td> an d for A — )• the single energy and the global so- 
lutions become identical. The optimum value for A is 
chosen from the condition 1 < X 2 /xltd < 1-05. The de- 
scribed two-step fitting procedure may be repeated sev- 
eral times by adjusting the energy dependence of the 
global solution, for example by changing the resonance 
parametrization to obtain a better agreement between 
the global and local solutions. 



IV. GENERATING PSEUDO DATA 

Measurements of the spin observables defined in Sec. 
|II C[ over a large angular range with accuracy sufficient 
to have an impact on multipole analyses, require exper- 
iments of high intensity, with linearly and circularly po- 
larized photon beams, strongly polarized targets and ac- 
curate recoil polarimctry together with hermetic detector 
systems. Such experiments have become technically pos- 
sible only very recently at the tagged photon facilities at 
ELSA (Univsersity of Bonn), CEBAF (Jefferson Labora- 
tory) and MAMI (University of Mainz). 

For this study, we have generated pseudo data for sin- 
gle and double polarization observables of the jp — > ir a p 
and 7p — > 7r + n reactions with statistical uncertain- 
ties comparable to those expected at the above men- 
tioned precision electromagnetic facilities within the next 
few years for energies from the reaction thresholds up 
to a photon beam energy of E 1 = 1500 MeV (W = 
1920 MeV) 



Events for the neutral and charged pion reactions are 
generated via a Monte Carlo algorithm. A sample of re- 
action parameters - the beam energy, meson scattering 
angles {9 and </>), circular and linear beam polarizations, 
longitudinal and transverse target polarizations and the 
recoil nucleon spin alignment angles (Or and 4>r) - are 
drawn from a weight function given by the polarized cross 
section, Eq.(14|-( 16). The events were generated for 
beam energy bins of AE 7 = 10 MeV and angular bins of 
Ae c,n _ 10Oj based Qn the jyjAID (MD07) model predic- 
tions. Beam polarizations are assumed to be Pt = 60% 
for linearly polarized photons and Pq — 70% for circu- 
larly polarized photons. These are typical values result- 
ing from e~ to 7 helicity transfer (giving circular photon 
polarization) and coherent bremsstrahlung (giving linear 
photon polarization), which are achieved at tagged pho- 
ton facilities. We assume that the beam polarizations 
have no energy dependence and therefore should be in- 
terpreted as resulting from various sets of experimental 
conditions. For target protons an average polarization of 
P = 80% is assumed, as is available with butanol frozen 
spin targets [16] . The polarization of recoiling nucleons 
may be measured using a subsequent scattering reaction 
in a carbon analyzer. For these measurements an aver- 
age analyzing power of A = 20% is assumed, close to val- 
ues typically achievable. For each observable, typically 
5 • 10 6 events have been generated over the full energy 
range. The pseudo data have not been folded with any 
additional acceptance given by a specific detector system. 

The pseudo data are an accurate representation of both 
the MAID MD07 ED solution and the observed reaction 
data, as shown in Fig. ([2]). Here, the unpolarized cross 
section, da/dQ, the photon beam asymmetry, E, the tar- 
get asymmetry, T, and the beam-target observable, G 
are shown at energies corresponding to the A(1232) res- 
onance. 

In the following section, we discuss the application of 
the SAID parametrization approach to obtain ED and 
SE solutions in fits to the pseudo data. We compare 
these solutions to the MAID MD07 ED solution and the 
pseudo data and discuss various features of their energy 
dependence. 



V. DEPENDENCE ON ENERGY 

Using the pseudo data generated from the MAID 
MD07 ED solution we apply the existing SAID 
parametrization approach 18, 19 to obtain the pion pho- 
toproduction multipoles. The results given below cover 
values of E 7 from 200 MeV to 1.2 GeV. The determina- 
tion of a set of isospin-1/2 and 3/2 multipoles in the low- 
energy region is problematic due to the distinct thresh- 
olds for ir°p and ir + n photoproduction from the proton. 
We avoid this issue by fitting data above the ir + n thresh- 
old. 

The use of Watson's theorem, greatly simplifies SE fits 
at energies where it is valid (depending on the multipole, 
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Multipole 


SE4 


MD07 


SP09 


Ml' 2 
M\' 2 
El' 2 
Ml' 2 
E 3 ' 2 
Mil 2 


9.40(0.08) 
18.06(0.16) 
2.28(0.15) 
9.26(0.27) 
1.82(0.05) 
2.52(0.18) 
1.00(0.07) 
51.91(0.06) 


9.36 
17.91 
2.21 
9.31 
1.79 
2.55 
1.12 
52.0 


7.30 
15.87 
1.65 
7.89 
1.76 
1.99 
1.08 
55.25 



TABLE II: Single-energy fit (SE4) to MD07 pseudo data compared to MD07 ED values at 340 MeV. The SP09 solution, fitted 
to the SAID database, is displayed for comparison. Multipoles given in millifermi units. 




FIG. 6: Multipole phases from energy-dependent fits. SP09 (solid), MD07 (dashed), ED4 (dotted). 



the Watson regime is typically somewhat higher than the 
first inelastic, two-pion threshold). Here, the multipole 
phase is fixed to the corresponding irN elastic partial 
wave phase and only the modulus allowed to vary. The 
SAID ED fit form, for each multipole, is 

M = a(l + iT„ N ) + PT^n + 7(IniT^ - \T^ N \ 2 ), (34) 

T-nN being the associated ttN partial-wave amplitude, 
which allows a smooth transition from the Watson 
regime. Here a contains the Born contributions plus a 
(real) phenomenological term, (3 (also real) is purely phe- 
nomenological, and 7 is a complex polynomial allowing a 
further departure from the Watson regime, proportional 
to the wN reaction cross section. 

As shown in Table |H| the SAID SE fit to MAID 



pseudo data, at the A (1232) resonance energy (_E 7 = 340 
MeV), reproduces the dominant amplitudes. The SE re- 
sult is, however, significantly different from an ED fit 
(SP09) 23J to the full SAID database of experimental ob- 
servables, covering the resonance region. We note that 
this database, determining the SP09 solution, includes 
few measured polarization observables. In producing the 
SE result, SE4, only da/dfl, P, S, and T have been fit- 
ted. 

In order to obtain an ED solution on which to base 
the SE solution of Table |Hj the SP09 solution was used 
as a starting point, from which a modified ED result was 
obtained from the fit to pseudo data, up to a photon 
energy of 1.2 GeV. A SE solution was then obtained, 
starting from ED, as described above. (Above, and in 
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Multipole 


nN PW 


MD07 


SP09 


ED4 


SE4p 


SE8p 




Sn 


16.7 


18.4 


17.4 


16.0(3.0) 


16.2(0.9) 


Ml' 2 


Pu 


72.7 


64.4 


73.2 


68.2(2.4) 


73.4(1.6) 


Ml' 2 




163.9 


167.1 


172.0 


176.8(6.5) 


167.5(1.8) 


El' 2 


£>13 


27.6 


17.9 


31.4 


27.5(3.1) 


26.1(1.0) 


mIL 2 


Dl3 


26.8 


22.2 


26.7 


25.5(1.8) 


26.7(1.0) 



TABLE III: Multipole phases (degrees) from single-energy fits to 4 (SE4p) and 8 (SE8p) observables at 600 MeV (phases 
searched), compared to the energy-dependent fits SP09, ED4, and MD07 (see text). Also listed is the associated nN partial 
wave (PW). 




FIG. 7: Real and imaginary parts of the Pn(pM) partial wave amplitude (Mj_ ). The solid (dashed) line shows the real 
(imaginary) part of the MD07 solution, used for the pseudo data generation, (a) Solid (filled) points display the SE4 fits to 
da/dQ and the three single-spin observables £, T and P. (b) Notation as in (a) for SE8 fits also including the beam-target 
spin observables E, F, G and H. Standard SE fits displayed; phases not searched. 



the following, the notation SEn or EDn denotes a fit to 
n observables.) We explore these solutions (the ED so- 
lution and its associated SE solutions) in detail below. 
The x-squared for the resulting SE fit (SE4) at 340 MeV 
is 116 for 144 dcr/dfl, P, S, and T pseudo data. We note 
here, for subsequent discussion, that only the I = 0, 1 
multipoles depart significantly from a Born approxima- 
tion at the A (1232) resonance energy. 

An alternative method that can be used at this en- 
ergy, which reduces the influence of Watson's theorem, is 
to first fit the multipoles for Tr + n photoproduction (both 
real and imaginary parts) for £ = 0, 1, using the real Born 
multipoles for higher waves, to fix the overall phase. One 
can then fit the 7r°p data separately, again fitting both 
the real and imaginary parts of contributing multipoles, 
with the Mf + P multipole phase fixed, given the associ- 
ated ttN phase (by Watson's theorem, but in this P33 
partial wave only) when combined with the previously de- 
termined MJi" multipole. This method requires fewer 
constraints, but results in multipoles with errors much 
larger than those given in Table |n] 

A form of this two-step method was used in an analysis 
by Grushin[20] and a recent update[2T] and was found to 
give multipoles with phases consistent with Watson's the- 
orem over the photon lab energy range of 280 < E 1 < 420 



MeV. As the P33 7rV partial wave remains essentially 
elastic for CM energies W < 1450 MeV (E 7 < 650 MeV) 
the method could, in principle, be extended to this en- 
ergy. This defines a second energy region, where not 
all multipoles have phases set by Watson's theorem, but 
reliable results can be obtained, with a self-consistent 
database of experimental reaction observables. A study 
by Omelaenko 22J suggests that the number of angular 
measurements required in order to determine multipoles 
up to a given angular momentum is of order 2£ for the 
differential cross section, and each single-polarization ob- 
servable: P, E, and T. An additional double polarization 
measurement, such as F or G, is required to fix a remain- 
ing discrete ambiguity. Here we are fitting 18 angular 
pseudo data for each observable, which clearly exceeds 
the requirements given in Ref. [22]. Note also that the 
case considered in Ref. [22] was neutral pion photopro- 
duction, and did not depend on the general use of Wat- 
son's theorem nor the charged-pion channel to determine 
the overall phase. 

In Fig. ([3]) , the SE and ED fits to pseudo data are com- 
pared to the cross section and single-spin asymmetries at 
340 MeV. Agreement of the SE and ED solutions with 
the pseudo data for da / dfl and the single-polarization ob- 
servables, P, S, and T is obtained; agreement of the pre- 
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FIG. 8: Real and imaginary parts of the Su(pE) partial wave amplitude (Eq+). Notation as in Fig. [7] 



dieted double-spin asymmetries (not shown) is obtained 
at a similar level. This has been accomplished without 
fitting the complete set of experiments required for am- 
plitude reconstruction, as described in Ref. [15) . 

In Fig.Q, this process has been repeated at 600 MeV. 
The SE fit to da/dQ, and single-spin asymmetries (P, S, 
T) at this energy agree with the pseudo data, with a x~ 
squared of 89 for 144 data. At this photon (CM) energy, 
600 MeV (1420 MeV), Watson's theorem is not generally 
valid in all multipoles. Inelasticity in the irN elastic P33 
partial wave is, however, small up to 1450 MeV[l]. If 

we assume the P 33 (pM) (M^| 2 ) multipole has the same 
phase as the ttN elastic P33 amplitude and that the high 
partial waves for 7r + n photoproduction are given by the 
Born terms, which are real, this should allow the recon- 
struction of all relevant multipoles. This was the ap- 
proach taken in Ref. [21] , although at lower energies than 
600 MeV. This expectation is supported by Figs. ^ and 
^ where the agreement between the displayed curves is 
of fair quality, though lower than that of the Fig.Q. 

In Fig. ([5]), for example, the prediction for beam-target 
observables from the SE fit is compared to the pseudo 
data and the ED solution (ED4). The agreement is gen- 
erally good but at the level of several standard deviations 
with respect to the pseudo data error bars at some an- 
gles. These discrepancies are reflected in the predicted 
multipoles 23J. 

The SAID SE procedure requires that the first ED fit 
to pseudo data provides the proper phases, as the subse- 
quent SE search holds these phases constant. We can im- 
prove the agreement, however, by relaxing the phase con- 
straint on the SE solution. In Fig. we compare phases 
of the initial solution (SP09) to the ED re- fit (ED4), and 
the solution underlying the pseudo data (MD07). The 
agreement between MD07 and ED4 is generally good up 
to about 600 MeV, though there are very large differences 
between MD07 and SP09 in some multipoles. 



In Table III the phases of several multipoles are com- 



pared at 600 MeV. Of particular interest are the multi- 

1/2 1/2 

poles E 0+ and E 2 _ , connected to the irN partial waves 
Sn and -D13, which are known to differ significantly be- 
tween MAID and SAID, as can be seen comparing the 
values from MD07 and SP09. Note that the ED fit has 
improved the agreement, but does not match the phases 
for the D13 multipoles, as shown in Fig. ([6]), panel (d). 
The SE fit to cross section and polarization observables, 
displayed in Figs.Q and ([5]), retains this phase mis- 
match. In the fit SE4p, the same set of four observables 
(da/dSl, P, S, T) is fitted, while allowing the multipole 
phases to vary; the result is tabulated in Table |III| In 
SE8p, the set of fitted pseudo data has been expanded 
to contain beam-target double polarization quantities (E, 
F, G, and H), again allowing the phase to vary. By allow- 
ing phases to vary, we can now obtain a good reconstruc- 
tion of the underlying model (MD07) for the dominant 
multipoles[23]. 



The 'standard' SE fit results (where the phases are 
not searched) are given for two multipoles in Fig.([7|) and 
The left and right panels display the effect of adding 
beam-target observables to the database. While uncer- 
tainties and scatter diminish with additional observables, 
some slight systematic deviations are visible. As de- 
scribed above, these can be reduced by allowing phase 
variation in the fit. 



At much higher energies, where multipoles connected 
to the nN P33 partial-wave are no longer constrained 
in phase, the agreement deteriorates. However, the 
phase comparisons in Fig. ^ suggest that the ED refit to 
pseudo data generated from MD07 obtains the approxi- 
mate phase of MD07 in many partial waves. It would 
be interesting to repeat this comparison with models 
other than MAID and SAID to better gauge the model- 
dependence of fit results at higher energies. 
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VI. CONCLUSION 

We have performed ED and SE fits to realistic pseudo 
data generated from the MD07 ED multipolc analysis 
over a range of photon (CM) energies from 200 (1120) 
to 1200 (1770) MeV. As a prelude, we have provided a 
detailed explanation of the amplitude definitions and the 
methods used in SE analysis by the MAID and SAID 
groups. This should be useful to those planning experi- 
ments and wishing to compare with the MAID and SAID 
predictions. It should also clarify the connections be- 
tween SE and ED fits, and the data they analyze. 

We have conducted a number of fits in order to ex- 
plore the model-dependence of SE results in specific en- 
ergy regions associated with pion photoproduction. Un- 
like the amplitude reconstruction process, which does not 
yield unique multipole solutions, the method described 
here appears to give good results with only fits to high- 
quality cross section and single-polarization data for en- 
ergies where Watson's theorem is valid. A second re- 



gion, extending over energies where the P33 ttN ampli- 
tude remains elastic, has also been found to yield repro- 
ducible multipoles. At higher energies, results become 
more model-dependent, and further study is required to 
determine the reliability of SE multipoles determined via 
the methods described above. 
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